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We propose a new type of Josephson junction formed by two superconductors close to the 
superconductor-Mott-insulator transition, one of which is doped with holes and the other is doped 
with electrons. A self-organized Mott-insulating depletion region is formed at the interface between 
two superconductors, giving rise to an asymmetric response of current to the external voltage. The 
collective excitations of the depletion region result in a novel phase dynamics that can be measured 
experimentally in the noise spectrum of the Josephson current. 



PACS numbers: 85.25.-j, 74.81. Fa 73.43.Nq 

The Josephson effect is one of the most fundamen- 
tal effects associated with the superconducting phase, 
regardless of differences among various superconduct- 
ing materials. In strongly correlated electron sys- 
tems, many other quantum phases exist in the vicin- 
ity of the superconducting phase, which gives rise to 
many salient features in Josephson junctions. For ex- 
ample, a long range proximity effect in superconduct- 
ing/antifcrromagnctic/superconducting (SAS) junctions 
was predicted [l| based on the competition between the 
superconducting and antiferromagnetic phases in high 
temperature superconductors @, 01, 0, [a]. An abrupt 
change of the Josephson critical current was predicted 
in the junction arrays near the superconductor-Mott in- 
sulating phase transitions. 

In this letter, we investigate a new design of the 
Josephson junction by taking advantage of the compe- 
tition between the superconducting and Mott-insulating 
phases. The two sides of the junction are hole and 
electron-doped superconductors, respectively, which are 
close to the superconductor-Mott-insulator transition. 
At the interface, a self-organized Mott-insulating region 
is formed as the tunneling barrier. We dub this junc- 
tion " Josephson diode!'' . Similarly to the conventional p- 
n junctions in semiconductor diodes, the depletion region 
is suppressed by a positive bias voltage and elongated by 
a negative bias voltage, giving rise to an asymmetric re- 
sponse of the Josephson current to the external voltage. 
While the geometry of the junction is similar to the p-n 
junctions in semiconductor diodes, there are fundamental 
differences. In the Josephson diode, the depletion region 
(Mott-insulator) is formed completely due to the quan- 
tum nature of competing phases. Moreover, the depletion 
region not only plays a role of tunneling barrier, but also 
its quantum fluctuation reveals important information of 
strong correlation effects. The fluctuation of the region 
boundaries couples to the carrier recombination process 
which results in an additional phase dynamics that can 
be measured experimentally in the noise spectrum of the 
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FIG. 1: The mean-field phase diagram of the Bose-Hubbard 
model. The left (right) side of the junction corresponds to 
the states at point A(B). 

Josephson current. 

To illustrate the physics of the Josephson diode, we use 
the Bose-Hubbard model for charged bosons. The long 
range interaction between bosons is approximated by the 
self-consistent effective potential V as 

(ij) i 

+ J2(V(i)-»)tf(i)b(i), (1) 

i 

where t is the hopping integral, U is the short range 
on-site repulsion. In homogeneous systems, V[i) equals 
zero due to the charge neutrality maintained by the back- 
ground charge, and the phase diagram is well-known as 
shown in Fig. [T] Mott-insulating phases appear at com- 
mensurate fillings and small values of t/U , which can be 
doped into superconducting phases either by particles or 
holes. 

The structure of the Josephson diode is depicted in 
Fig. [5] a: its left and right hand sides correspond to two 
points A and B in Fig. [TJ where the particle densities 
(n(i))o — 1 =F 8 for i on the left (right) side, respectively. 
Due to the chemical potential imbalance, bosons diffuse 
from the right to left side. Since a Mott-insulating re- 
gion interpolates between the hole-doped superconductor 
A and the particle-doped superconductor B, it appears 
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FIG. 2: The sketch of the junction: (a) the spatial distribution 
of chemical potential before equilibrium is formed; (b) the 
Mott-insulating depletion region and charge distribution in 
equilibrium; (c) the potentials Vh and Vd for the holons and 
doublons across the junction, respectively. 



across the junction in the real space as the depletion re- 
gion depicted in Fig. [2]b. We define /ii(2) as the chemical 
potentials at the boundary of the Mott-insulating phase 
depicted in Fig. [TJ Inside this region, the density (n(i)) 
is fixed to the commensurate value 1 and the charge neu- 
trality is no longer kept. The resulting internal electric 
field E gives rise a change of the electric potential AV 
across the junction to compensate the local chemical po- 
tential difference. As a result, the total chemical poten- 
tial jj, becomes constant. The E field is given by 



WE(i) = -VY(i) = |(n(i) 



(2) 



where q is the charge of bosons, e is the dielectric constant 
and the lattice constant is taken as 1. In the following, we 
consider the case that A/x = — Ml is much larger than 
/til — fi'A and [is — (J>2 ■ In this case, most of the charge non- 
neutral region is the Mott-insulating, and we neglect the 
contribution to the E field outside the Mott-insulating 
region. Then Ay across the depletion region can be 
approximated as AV^ « — J_ D E{x)dx — D qS = A/i, 
where 2D is the length of the Mott-insulating depletion 
region, thus 
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Next we discuss the Josephson diode effect under an 
external voltage V ex . In the conventional Josephson junc- 
tion, the coupling term is given by 



- Jcosf 
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(4) 



where the coupling strength J is not sensitive to V ex . 
However, the situation in the Josephson diode junction 
is quite different. When V ex is applied, it not only cre- 
ates the time dependent phase difference across the junc- 
tion, but also changes the length of the depletion region. 
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FIG. 3: The length of Mott-insulating depletion region D 
and the Josephson coupling strength J as a function of the 
external voltage V ex ( arbitrary unit for the vertical axis). 



Similarly to the case of semiconductor diodes, the de- 
pletion region is suppressed (elongated) in the forward 
(reversed)-biased junction since the external electric field 
aligns in the opposite (same) direction of the internal 
electric field. By a similar reasoning for Eq. [3l we obtain 
D in the presence of V ex as D m ^/e(A/i — V ex )/(q5). 
When D is much longer than the coherent length £, the 
Josephson coupling J is proportional to e~ a ~ , where a 
is a dimensionless constant. Therefore, J as a function 
of V ex is asymmetric as 



log J(V ex ) cc yj (A/i - V ex )/(qS). 



(5) 



which is sketched in Fig. [3] This is a natural generaliza- 
tion of the effect of the p-n junction in semiconductors. 
Please note that when V ex is close to A/i, the Mott de- 
pletion region is almost completely suppressed, then Eq. 
[5] does not apply. 

A major difference between the Josephson and con- 
ventional semiconductor diodes is that the charge cur- 
rent is no longer unidirectional in the former case. In the 
semiconductor diode, the unidirectional charge transport 
clearly breaks time reversal symmetry, which is consistent 
with its dissipative nature. In contrast, in the Josephson 
diode, the dissipationless transport keeps time reversal 
symmetry. As a result, in the zero bias case, the DC 
Josephson current can flow either from the p to n-side or 
equivalently well from the n to p-side. At the finite bias, 
the AC Josephson current flows back and forth between 
the two sides. The asymmetric response to the external 
voltage lies in the magnitude of the Josephson critical 
current, which is an effect from the explicit parity break- 
ing in the structure of the Josephson diode junction. 

In a semiconductor p-n junction, there are the pro- 
cesses of the diffusion and recombination of charge car- 
riers. The analog of such processes also exists in the 
Josephson diode junction. We introduce the singlon 
(st(i)), holon (ajj(i)), and doublon (aj,(i)) operators for 
the empty, single and double-occupied states on the site 
i, respectively. We limit the Hilbert space on each site 
to these three states, which is a good assumption for the 
Bose-Hubbard model when the density is close to one. 
In this case, the holons and doublons are considered as 
the majority carriers in the left and right sides, respec- 
tively. Since each site can only have one such boson, these 
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bosons are hard core-like, thus namely they satisfy the 
following constraint a h ah+s^ s+a p a p = 1. We can rewrite 
the original boson operators in terms of these hard core 
bosons as W — s^a^ + aps. The density operator of the 

original bosons is given by n = b^b — a p a p — a^dh + 1- 
The hopping term in this representation becomes 

H t = -*5^{[o^(j)ofc(i) + 4(*)ap0')]* t 0>(*) + ^c. 

+ K(j)a p (j) + a p {i)a h {j)]s\i)s\j) + h.c.X, (6) 

where the first two terms describe the hopping processes 
of holons and doublons and the second two terms describe 
their recombination and regeneration processes. Corre- 
spondingly, the current operator J cur is given by 

Jcur = it^2{[al(i)a h (j) - aj(j')a p (i)]s t (i)s(j) - h.c.} 

-it^2{[a h (i)a p (j) - a h (j)a p (i))s* (j)s^ (i) - h.c.}, (7) 

ij 

where the first two terms are the drift currents and the 
last two terms are recombination and regeneration cur- 
rent. 

In order to study the Josephson coupling explicitly in 
the diode, we take the following general mean field trial 
wavefunction, 

* = Yl {cos e^s^i) + e^ 2 ^ sm9 1 (i)[e lMt) 

i 

cos0 2 (i)4(i) + e-^M sin0 2 (i)4(i)]}, (8) 

which gives the superconducting order parameter as 

(6 t W) = sin2 ^(*) e ^(i)( e ^W cos g 2 (i) 

+ e-^ 2 «sin6l 2 (i)). (9) 

There are two independent phase variables 0i j2 in the 
above formulae. </>i(i) is the conventional phase of su- 
perconducting order parameter, which is the conjugate 
variable to the density, namely, [n(i), </>i(i)] = i- The 
other variable, </> 2 (i) is associated with the quantity 
a p a p + a^ah = 1 — s^s, namely, the density of singlons 
n s (i) = s's. Thus, we have [n s (i), 4>2{i)\ = i- In the con- 
ventional Josephson junction, <p2{i) can be simply fixed 
to be zero since the number of singlons can arbitrarily 
fluctuate. However, the presence of the depletion insu- 
lating region changes the above picture. Not only does 
the depletion region serve as a tunneling barrier, but also 
its fluctuation leads to additional physics, which is the 
fluctuation of the number of singlons. Assuming the vari- 
ables of the wavefunctions in the superconducting states 
in the left and right sides of the junctions are given by 



{(t>u,r, <p2i,r, 8il,r, 02i,r), we rewrite Eq. [6] as 

H t = Jicos(A0i - A0 2 ) + J 2 cos(A0i + A0 2 ) 

+ J 3 cos(A0i -tp) + J 4 cos(A</>i + tp), (10) 

A(f>i = A(j> ir - A(j) a (i = 1, 2), tp = fci + 2r , 

Ji = sin 2#i/ sin 2Q\ r cos 621 cos # 2r / 4, 

J 2 = sin 2^1/ sin 26*i r sin0 2 ; sin0 2r /4, 

J3 = sin 2#i/ sin 2Q\ r cos 821 sin 2t -/4, 

J4 = sin 28u sin 29\ r sin^ cos# 2r /4. (11) 

Similarly to the Josephson coupling in multi-band su- 
perconductors 0, B @| , Eq. [10] indicates hidden internal 
dynamics in the junction. There are two observations. 
The first is that the coupling depends not only on A</> 2 , 
but also on tp, the sum of the phases <p2i and 2r , reflect- 
ing the recombination and regeneration processes. As 
shown in Fig. [TJ holons and doublons are the major- 
ity carriers in the left and right sides respectively, thus 
cos6* 2 / ^> sin6> 2 ; and sin# 2r 3> cos# 2r hold. As a result, 
the J3 term describing the recombination and regenera- 
tion processes dominates over other terms. The second 
is that if cf>2 is set to 0, we recover the conventional form 
of the Josephson coupling energy i.e. H = Jcos(A^i). 
Without the insulating region, the number of singlons is 
not determined in both sides. Therefore, tp and A</> 2 do 
not have any dynamics and can be chosen to 0. How- 
ever, their dynamics can be created in the presence of 
the depletion insulating region. 

We discuss the fluctuation of the Mott-insulating de- 
pletion region. When the system is slightly out of equi- 
librium, the depletion region has to shrink or expand 
symmetrically respect to the center of the junction in or- 
der to maintain the charge neutrality. This corresponds 
to the recombination and regeneration processes in the 
depletion region, which changes the value of nf + nf, 
and thus creates the dynamics of tp. The frequency 
of this mode can be estimated as follows. The mis- 
match of chemical potentials in the two sides of the junc- 
tions depends on the length of the depletion region as 
Afjbmis = A^i — Ay = Afi - D 2 qd/e. The restora- 
tion force on the charge carriers can be approximated 
as F = —Afi mis /(2D) = -AxAfi/D 2 , where Ax is 
the small displacement of the boundaries of the deple- 
tion region. This mode is essential the plasmon mode 
with the frequency determined as to 2 = qd/(m*e), where 
to* sa l/{tS) is the effective mass of the charge carrier. 
Thus tub scales linearly with the doping level d. 

In general, the fluctuations of the phase tp is a diffu- 
sion process. Let T be the decay rate set by the diffusion, 
which depends on disorder and detailed low energy ex- 
citations in the insulating and superconducting regions. 
The fluctuation of tp lead to the following correlation 
function, 

/ e i<p(,i) e -iy(o)\ _ e iuj b t-r\t\ 
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Therefore, the total current correlation function , 

(J(t)J(O)) oc cos(a;e,t)e~ r|t| , (13) 
which gives the noise spectrum of current as 



■S-UM = j dte^{{J(t)-(J(t)),J(0)-(J(0))}) 

r r 



(14) 



Now we briefly discuss the case for neutral bosonic sys- 
tem, say, the p-n junction made by cold bosonic atoms. 
The Bose-Hubbard model [1(3] has already been realized 
in optical lattices experimentally. In such systems, it re- 
quires an external potential drop to create a depletion re- 
gion in an inhomogeneous optical lattice. In this neutral 
system, there is an additional mode, a sliding mode, as- 
sociated to the collective excitations of depletion region. 
For the sliding mode, we can imagine that the system is 
slightly out of equilibrium by sliding the whole depletion 
region. Since this mode changes the value (nf — rig), it 
naturally creates the dynamics of A</>2. In the two band 
picture, this mode is corresponding to transfering a hole 
(or doubleon) in the one side to the other side. Since 
the difference of band energy is given by the mott gap, 
Afi which is proportional to U, the mode is expected to 
have the frequency uj a = A/i. This mode has the same 
response in the noise spectrum of the junction as the 
plasma mode expect their energy difference. A design of 
the "atomtronic diode" exhibiting asymmetric response 
to external chemical potential difference was discussed by 
Seaman et al. [ll|. The major difference to our design is 
that no insulating depletion region exists in their design. 

Experimentally, Bose-Hubbard mode for charged 
bosons can be realized in Josephson arrays By 
choosing the chemical potential configuration discussed 
above, the Josephson diode should be formed. High T c 
cuprates are another class of possible systems to realized 
the Josephson diode. Two different insulating phases ex- 
ist in cuprates. One is in undoped systems and the other 
is at the 1/8 doping level [12, EH, both of which can 
be used to create the Josephson diode in principle. A 
universal asymmetric density between positive and neg- 
ative bias voltage in STM experiments 0, EE El has 
been observed, which can been explained by considering 



an insulating competing order at surface [Fjj, IS, 3 2Cj |- 
This is an important support for the possible realiza- 
tion of the proposed junction, although other alterna- 
tive explanations 2l|, |22j also exist. The key to cre- 



ate the Josephson diode junction is a spatially selective 
doping control in cuprates. An in-plane SAS junction 
has been created by using the spatially selective and re- 
versible doping control techniques in cuprate films [23[. 
Photodoping tec hniq ues can also locally control the dop- 
ing level [3, HE l26l |. These techniques can be used to 



create two close-by regions with different doping levels 
with an insulating phase between them. Although it is 
clear that the Josephson diode in cuprates can not be 
described by the simple Bose-Hubbard model, we believe 
that the physics presented here is still valid. The effec- 
tive bosonic Hamiltonian in cup rates has to include the 
spin degree of freedom [3, [E [27( . Since we are interested 
in the charge transport in the Josephson diode, the spin 
degree of freedom does not have direct effect on the prop- 
erties we studied in this paper. However, it is largely an 
open question whether any other new properties can be 
induced by the spin degree of freedom. 
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